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Abstract. Positional ordering of a two-dimensional fluid of hard disks is examined 
in such narrow tubes where only the nearest-neighbor interactions take place. Using 
the exact transfer-matrix method the transverse and longitudinal pressure components 
and the correlation function are determined numerically. Fluid-solid phase transition 
does not occur even in the widest tube, where the method just loses its exactness, but 
the appearance of the dramatic change in the equation of state and the longitudinal 
correlation function shows that the system undergoes a structural change from a fluid 
to a solid-like order. The pressure components show that the collisions are dominantly 
longitudinal at low densities, while they are transverse in the vicinity of close packing 
density. The transverse correlation function shows that the size of solid-like domains 
grows exponentially with increasing pressure and the correlation length diverges at 
close packing. It is managed to find an analytically solvable model by expanding the 
contact distance up to first order. The approximate model, which corresponds to the 
system of hard parallel rhombuses, behaves very similarly to the system of hard disks. 

Keywords: confined systems, positional order, transfer matrix method, exactly solvable 
models 
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1. Introduction 

The knowledge of the thermodynamics of fluids in extreme confinements has a great 
importance due to the wide range of chemical and biological applications and the 
appearance of nanoscale-fluids in nature [I]. The most important examples are the 
micro- and nanofluidic devices (2j [3], carbon nanotubes [U [5], biological ion-channels 
[6] and adsorption in porous materials such as the zeolite [7]. To capture the essence 
of the rather complex behaviour of such confined systems it is important to apply some 
simplifying assumptions for the particle-particle and particle-confining wall interactions. 
Using hard body molecular models and hard confining walls it is possible to study purely 
the geometrical packing effects as the energetic effect are not included in these models. 
In this regard much attention has been paid to understand the structure and the phase 
behaviour of several confined systems such as the ordering of hard spheres in narrow 
cylindrical pore [8j [9j [10] and slit pores [11] and the helical packing of soft spheres in 
nanotubes [T2| IT3] . The rapid development of the X-ray techniques makes also possible 
to study extremely confined fluids, where the particles are confined between two hard 
walls with thickness just slightly higher than the diameter of the particle [14J. These 
systems can be often modeled as quasi-one-dimensional systems. 

Since the seminal work of Tonks [15 J considerable progress has been made to 
understand the thermodynamic and structural properties of one- dimensional systems. 
The pair correlation function was determined in various systems interacting with short- 
range pair potentials [T6l [T71 HE], the link between the short-ranged pair potential and 
the pair correlation functions was examined in a charged- stabilized colloidal system 
trapped to one dimension [19] and the long-range interaction induced freezing transition 
of hard rods was studied by [20]. Even our knowledge of the ordering properties of one- 
dimensional hard rod fluid (Tonks gas) has been widened by the exact determination of 
the percolation size [2T] and by the study of the residual multiparticle entropy at the 
onset of solidlike behaviour [22J. 

The examination of confined two-dimensional (2D) hard particle fluids can be 
considered as a natural extension of the one- dimensional studies, if the available space 
for the particles is infinitely long in one dimension, but it is very restricted in other 
dimension. This condition guarantees that the phase behaviour of the confined fluid 
does not deviate substantially from that of Tonks gas. Of course, with weakening 
confinement the system gets naturally closer to the two-dimensional bulk fluid and the 
thermodynamic properties of such fluids get very far from that of Tonks gas [23, 24J. 
Several studies are devoted to extend the exact Tonks theory of one-dimensional fluid 
to the confined two-dimensional one [25, 26, [2TJ, EE]. If the confinement and the pair 
interaction does allow only nearest neighbor interactions, it is possible to derive such an 
integral equation which provides the exact thermodynamics of the system. In the system 
of 2D hard disks confined into a narrow tube, the equation of state has been determined 
by the iterative solution of the corresponding eigenvalue equation [27]. Later analytical 
equation has been derived for the equation of state in the high pressure limit, which is not 
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Figure 1: Schematic representation of the system of hard disks confined between two 
hard walls, where D is the diameter of the disks and W+D is the wall-to-wall separation. 
Overlaps are not allowed between the disks and the vertical position of the disks is 
confined in the interval of — W/2 < y < W/2. 

valid at low densities [29|I30|. Recently both the transverse and the longitudinal pressure 
components of the hard and soft disk have been determined at such low densities, where 
the second virial theory is still valid [3T| [32] . In this paper we extend the previous 
studies to examine the local structure of the hard disk fluid at both low and high 
pressures and we devise an analytically solvable model system which follows closely 
the phase behaviour of the hard disks. We calculate the transverse density profile, the 
pair correlation function and the longitudinal pair distribution function to examine the 
growth of the solid-like clusters with increasing pressure. We pay special attention to 
the effects of the structural change from a Tonks gas behaviour to a solid-like one. 

Our paper is organized as follows: the transfer matrix method is presented for 
the confined hard disks in the next section. The eigenvalue equation is derived for the 
determination of the Gibbs free energy, the equation of state, the correlation function 
and the local density. Moreover we present the details of our MC simulation study 
for the pair correlation function measured in the direction of the tube' s long axis. In 
section [3] we show that analytical results can be derived for both thermodynamic and 
structural properties if the contact distance is expanded up to the first order. The 
results are shown to be exact for the system of parallel hard rhombuses placed into 
a very narrow tube, but it is approximate for hard disks even in the vicinity of the 
close packing limit. We compare our numerical, analytical and MC results and we also 
examine the fluid-solid like structural change, which takes place with increasing density 
in section HI Conclusions are given in section [5j 

2. Transfer matrix method of the confined two-dimensional hard disks 

We consider an isobaric ensemble of N freely moving two-dimensional hard disks in such 
a restricted geometry where only nearest neighbor interactions are allowed (see figure 
[1]). The interaction potential between two neighboring hard disks is defined by 




(1) 
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where r iti+ i = \r i+ i — is the length of the centre-to-centre distance vector between 
particles i and i + 1, = (xi,yi), Xi and x/i denote cartesian coordinates of the position 
vector of the particle i and D is the diameter of the disk. Equation (CQ) precludes the 
overlap between the particles. To restrict the number of interactions only to the nearest 
neighbors a confining external field is applied, which acts only on the y coordinates of 
the particle's centers and allows the particles to stay in a very narrow tube such that 

( x / 0, \ Vi \ < W/2 
UUm) = { oo, \ Vi \ > W/2 ' (2) 

where W is practically the width of the tube even if there is an additional D distance 
in the tube which is inaccessible for the centers of the particles (see figure 1). It is easy 
to show that the tube width must be in the interval < W < D\f?>/2 to exclude the 
next-nearest neighbor interactions. In such a narrow tube the system is anisotropic, 
i.e. the horizontal and vertical pressure components P xx and P yy are different. From 
the two possible options to define the Gibbs free energy (G?) from the Helmholtz free 
energy (F) we choose the horizontal pressure component to get G from the Legendre 
transformation, i.e. G = F + P XX A, where A = LW is the accessible area for the centre 
of the particle and L is the length of the system along the x axis. The corresponding 
isobaric partition function of the confined system can be written as 

i 7 7 ( N \ 

z np xx t = J dxi ■ ■ ■ dx N J dyi--- dy N exp I -/3 {u(r i>i+ i) + P xx Wx iti+1 ) J , (3) 
o -W/2 ^ i=1 ' 

where = l/k^T (T being the temperature and ks the Boltzmann constant) and 
A = hj \j2-nk^Tm is the thermal de Broglie wavelength. Note that the missing iV! term 
in the partition function is due to the use of the fixed x% < x<i < ■ ■ ■ < x^ order. This 
order can be fixed because the tube is enough narrow. Furthermore the finite integral in 
Xi is replaced by the integral with infinite upper bound in the isobaric ensemble and we 
use the fact that L = Y^Li where x^j+i = Xi+\ — X{. Changing from the integration 
variable Xi to x^+i and using equation (JTJ) it is easy to prove that 

W/2 



z np xx t = J dyi ■ ■ ■ dy N 

-W/2 



J dxi )2 exp (-(3P xx Wx 1<2 ) ■■■ J dx Ni i exp {-/3P xx Wx N>1 ) , (4) 



where o~(yi,yi+i) = J D 2 — (jji — yi+i) 2 is the horizontal projection of the contact 
distance between two hard disks being at the positions yi and yi+i. The lower bound of 
the integrals in Xi comes from the overlap condition, which prescribes that r^+i > D. 
We have exploited the periodic boundary condition in the x direction, too, which means 
that r^v+i = ri- Now we are in the position to perform the positional integrations in 
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X{j + i analytically, i.e. 

TV W / 2 



J NP XX T 



(a 2 , 



) J dj/i • • ■dy N K(yi,y 2 )K(y 2 ,y 3 ) • • • K{y N ,y x ) , (5) 

where K(yi, yi+i) = exp(—(3P xx Wa(yi, yi + \))/ (3P XX WD can be considered as an element 
of an infinite-dimensional matrix (K) . In this formalism the matrix product is given by 

W/2 

K 2 (y h y i+2 ) = J dy i+1 K(yi,y i+1 )K(y i+1 ,y i+2 ). Therefore the partition function can 

-W/2 

be considered as the trace of the matrix K N , that is Z^p^t — TrK N = J dyK N (y, y). 

-W/2 

Due to the fact that the result of trace operation is independent of the used basis, it is 
most convenient to evaluate the partition function in the eigenfunction basis where the 
matrix K is diagonal. Denoting the eigenvalues of K in decreasing order in the absolute 
value as |A | > |Ai| > |A 2 | . . ., equation (jSJ) reduces to Z NPxxT = X)£L Aj^, where the 
eigenvalues are the solutions of the following eigenvalue equations 

W/2 

dy'K(y,y')4> k (y') = X k i>k(y) , (6) 

-W/2 

where ipk{y) is the kth eigenfunction (k — 0, 1, 2, . . .) and K(y, y') is the integral kernel 
defined in equation (j3J). To our best knowledge, the above eigenvalue equation does not 
have analytical solutions for hard disks. Therefore we solve it by the iteration using the 
trapezoidal quadrature for the numerical integration. We can obtain the working form 
of equation fl6]) by the use of the normalization constraint of the eigenfunctions, that is 

W/2 

J dyipk{y) 2 — 1> In this way we write the equations of the eigenfunctions and the 

-W/2 

corresponding eigenvalues as follows 

W/2 

J 2 dy'K(y,y')My') 

My) = w/2 ~ w/2 , (7) 



/ dy'dy''K(y',y")My")My') 

-W/2 



and 



W/2 

X k = J dydy'K(y,y')My')My) ■ (8) 

-W/2 

In the thermodynamic limit (N — > oo) the largest eigenvalue (Aq) has the most 
significant contribution to the Gibbs partition function, i.e. the Gibbs free energy 
takes the following form: (3G/N = — ln(Ao/-D), where the constant term proportional 
to ln(A/D) is omitted. The equation of state can be obtained form the standard 
thermodynamic relation between the Gibbs free energy and the pressure, which is now 
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W/p = dg/dP xx , where p = N/L is the linear density and g = G/N. Similarly the 
transverse pressure can be also derived from the Gibbs free energy as follows 

dg 



Pyy Pxx P 



dW 



(9) 

P 

1 xx 



9Kh3. 




Regarding the positional distributions of the disks, ipo(y) eigenfunction gives 
account of the positional ordering along the vertical direction, because ipoiu) * s the 
probability density distribution function and p(y) = pip^y) is the particle number 
density distribution function along y axis. Other important function which measures the 
posititional correlation in the vertical directions between two particles % and j is defined 
as g(i,j) = (iyi — (yi})(yj — (%•})/• For the present case we only need to calculate the 
NP XX T ensemble average of the products of the horizontal positions of particles % and j 
(g(hj) = (yiVj))) because the average position of the particles is trivially zero ((t/j) = 0). 
It can be shown by using the eigenvalues and eigenfunctions that 

W/2 1 2 

J d v Mv)vMv) ■ ( 10 ) 

-W/2 

Note that the fluctuation in the vertical position is given by 

W/2 

9(i,i)= J dyij 2 (y)y 2 . (11) 

-W/2 

It can be seen from equation f fTUj) that the correlation function decays exponentially for 
\i — j\ — > oo, i.e. we can write that g(i,j) ~ exp(— \i — where the correlation 

length (£) is given by ^ = (In | Aq / A i | ) — 1 . Besides the transverse correlations it is also 
necessary to get information about the longitudinal positional order to characterize the 
structure of the confined disks. We do this by the MC simulation of the longitudinal pair 
distribution function. This function is computed by our own Monte Carlo simulation in 
canonical ensemble (NAT) from the expression 

g[x = \x - x"\ J = p~ 2 ( ^E^ 1 ' ~ x ') $( x j ~ x ") ) ■ ( 12 ) 

\ * j¥* ' NAT 

In this expression the configurations are collected from the simulation of 1000 hard 
disks, which are sampled from 10 5 configurations beyond an equilibration process of 10 5 
configurations. The starting configuration is a triangular lattice, where the particles are 
at the two walls in alternating up-down sequence. In addition, the periodic boundary 
condition is used in the horizontal direction in all studied cases. Using the Metropolis 
sampling, the particles are chosen and moved randomly with an adjusted maximum 
displacement to ensure the requirement of 40-50% acceptance probability of all trials. 
Due to the simple form of the pair potential and the external field, one trial move of a 
chosen particle is accepted if the particle does not overlap with its neighbors and the 
walls. 
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3. Analytically solvable model 

In this section we proceed by searching for the analytical solution of the eigenvalue 
equation ([6]). We are facing an integral equation which corresponds to a homogeneous 
Fredholm equation of the second type. Due the non-separable form of the contact 
distance o~(yi,yi+i) and the integral kernel, it is well-known that the Fredholm equation 
does not have a general solution. In previous studies, Kofke and Post [27J solved equation 
(EJ) by the standard iterative solution, while Kamenetskiy et al. [30] considered only the 
high pressure limit and performed a perturbative calculation. 

Here we show that equation ([6]) has analytical solutions in a closed form if the 
horizontal projection of the contact distance is given in the following form 

Vi+i) = c + c\vi - y i+1 \ , (13) 

where Cq and c are constants. This formula can be seen as an approximative expression 
for the contact distance of confined hard disks. Expanding the exact contact distance 
of disks around yi — yi + i = ±W up to first order in yi — y i+ i we obtain 

C ° ~ VD 2 - W 2 ' ° ~ VD 2 - W 2 ' [ ' 

The reliability of the first order Taylor expansion is presented in figure [2j where the 
curve of the exact contact distance is compared to the approximate formula ( [TBI at 
different widths of the tube. We can see that the equation (fT3|) works quite accurately 
for narrow tubes, but the error of ( fTBT) can even exceed 60% for W = 0.8D at zero 
distance. However, at close packing limit this expansion is justified even at wider tubes 
because the neighboring disks are forced to be at the opposite walls to fill in the space 
most efficiently Furthermore, equation (fTBl gives exact formula for the hard rhombuses 
which can not rotate and one of their diagonal axes is constrained to be parallel to the 
wall. Denoting the lengths of the transverse and longitudinal axes of a rhombus as a and 
b respectively, it is easy to show that a = b — b\yi — yi+i\/a. Now we are at the stage to 
present the solutions of equation ([6]), which is now supplemented by equation (fTBl with 
Cq and c < constants. The largest eigenvalue and the corresponding eigenfunction are 
given by 

2c / CqP* 

A = D . —-: o eX P FT - )' ( 15ft ) 




where k,q is the unique positive solution of the following equation 

( W \ 

kq = —cP* coth k — - , (15c) 



2D) 

and P* = (3P XX WD is the dimensionless preassure. Using the above eigenvalue we can 
determine the Gibbs free energy, linear density and the transverse pressure component, 
while the transverse density profile can be obtained with the aid of ipo eigenfunction. 
It is important to know the other solutions of the eigenvalue equation, because the 
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Figure 2: Horizontal projection of the contact distance between two disks as a function 
of their vertical distance. The first order approximation of the excluded distance is 
shown for different tube widths. 



pair correlation function (llOl) contains them. In absolute value the second largest 
eigenvalue, Ai, and the corresponding eigenfunction, ipi, depend on the pressure. In 
the high pressure case, when —cP*W/D > 2, 

A i = D 7 r,l? — 2 ex P ( - -tt) ' ( 16a ) 



(cP*y - k\ 1 \ D 

1 /2 

ih 1 (y) = -$ = ( ^ — ] sinhf/^i^-V (166) 

where K\ is the unique positive solution of the following equation 

( W \ 

K\ = —cP* tanh I k>i—— I . (16c) 



2D) 

In the low pressure case, when — cP*W/D < 2 
x 2c / c P* , 



{cP*f + «? ^ V D 

^^Tb Uw/d-Z^w/d) ) sin H)' (m) 

where K\ is the unique solution of the following equation in the interval K\ G [0, ttD/W] 
«! = -cP* tan f kj — J . (17c) 
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In the case of — cP*W/D = 2 the above two groups of formulas (16) and (17) give 
the same limit: X t = D2cexp (-c P*/D) /(cP*) 2 , and ipi{y) = 2V3W~ 3/2 y. We must 
mention that Ai < in all cases. This is the reason why the correlation function shows 
alternating behaviour. 

All of the remaining eigenvalues and eigenfunctions can be summarized for odd and 
even indices as follows. 



Xi = D 



2c 



My) 



( C P*) 2 + K 2 

2 f\}A 



cxp 



c P* 
D 

1/2 



D {KiW/D + siniKiW/D) 

2 k,: 



COS K 



D 



1/2 



D \K i W/D-sm(K i W/D) j 



sin Kj — 
D 



2,3,... 

if i is even 
if i is odd 



:i8a) 



;i86) 



where k^ (i — 2, 3, . . .) are the unique solutions of the following equation in the intervals 
Ki G [(i - l)irD/W,inD/W] 



K; 



w 

+CP* COt ( Ki— — 

2D 
W 



-cP* tan ( Kj 



'2D 



if i is even 
if i is odd. 



;i8 C ) 



Finally we mention that equation ([6]) has also analitic solution for c > 0, but they 
are different and have relevance for the system of freely rotating hard needles in one 
dimension 1 33 1 . 



4. Results and discussions 



In this section we begin with the presentation of our numerical results for the pressure 
components, density profiles and some structural properties. To do this we make a 
connection between the present model and an analytically solvable reference system. 
In the limit of W — > 0, the system becomes a one-dimensional fluid of hard rods with 
diameter D, which is the well- know Tonks gas [15J. We can achieve the Tonks gas 
behaviour at finite width, too, by equating the contact distance to D at any transverse 
separation, i.e. a(yi,yi + i) = D for all horizontal positions. This case corresponds 
to the system of parallel hard squares in a tube. It is easy to show that equation 
([6]) gives ipQ = 1/yW and Ao = exp(—/3P xx W D)/(/3 P XX D) for the squares and the 
relation between the pressure and the Gibbs free energy gives back the Tonks equation 
for the pressure, which is j3P xx W = p/(l — pD). From the equation of transverse 
pressure (Q we get that the squares behaves as an ideal gas between the two walls, 
i.e. PPyy = N/(LW). Since the contact distance deviates only slightly from D in 
very narrow tubes (see figure [2]), the confined parallel hard squares can considered as 
a reference system of the hard disks. Our numerical results show that the system of 
hard disks follows the phase behaviour of hard squares in narrow tubes. This can be 
observed even in figure [3] for W = D/2, where the transverse and longitudinal pressure 
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components are shown as a function of linear packing fraction (77 = pv ' D 2 — W 2 ). 
The longitudinal pressure follows the Tonks equation, while the transverse one can 
be described by the ideal gas law up to 80% occupation of the longitudinal direction 
(77 Rj 0.8). The transverse density profile also shows that the density is almost constant 
across the tube for /3P XX WD = 5, i.e. i/jq = l/y/W is a reasonable approximation up 
to i] w 0.8 for hard disks. For higher pressures (packing fractions) the density profile is 
peaked at the walls and there are less particles in the centre of the tube, which is shown 
at (3P XX WD = 25. The results are more interesting at the possible widest tube width 
(W = Dy/3/2), which can be studied exactly taking into account the nearest neighbour 
interactions only. The Tonks gas and the ideal gas behaviour are lost at much lower 
packing and the derivative of the longitudinal pressure shows non-monotonic behaviour 
with the density. It has two inflection points at intermediate packing fractions and it 
is crossed by the transverse pressure in the vicinity of larger inflection point. These 
results suggest that some positional rearrangement must take place in the system to 
give rise a big change in the pressure. Especially the cross point between the two 
pressure components is very informative. It means that the particles must interact 
with each other very strongly along the transverse direction, which can only happen 
if the particles are in alternating arrangement at the walls. This speculation is also 
supported by the density profiles shown in the inset of figure [3j A deeper understanding 
of the positional ordering can be gained from the transverse pair correlation function 
ffTU]) . which measures the propagation of the fluctuations between the particles. Figure 
H] highlights two important features of the positional ordering. Firstly the correlation 
length diverges with the pressure (density), which shows that the size of the ordered 
clusters increases and a solid-like structure may emerge at the close packing. Secondly 
the oscillatory behaviour of the correlation function gives information about the spatial 
arrangement of the disks, namely if one particle moves into one given transverse direction 
its neighbor does it in the opposite way. This type of correlation can occur only if the 
particles are located in an alternating up-down sequence. At low pressures only short 
range alternating arrangement takes place, but the ordered clusters grows exponentially 
at higher pressures and the arrangement become similar to a section of a triangular 
lattice. As the correlation length does not diverge below the close packing, no phase 
transition between fluid and solid phases can occur in the physical meaningful density 
range. We can also see that the transverse correlation becomes stronger in wider tubes, 
as the alternating packing becomes entropically more favored and it is getting harder 
to make jumps between the two walls. 

The longitudinal pair distribution function also helps us to understand the 
structural changes occurring in the tube. The results of our canonical MC simulation 
shows that in case of low pressures the average distance between the neighboring disks 
is in the order of D for both narrow and wide tubes. This means that at low packing the 
confined hard disk systems behave like the Tonks gas. Figure [5] shows that one additional 
first peak emerges in the pair distribution function at intermediate packing, which means 
that the system is just between the Tonks gas and solid-like behaviours. In the very dense 
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Figure 3: Horizontal (P xx ) and vertical pressure (P yy ) components as a function 
of packing fraction for W = D/2 (upper panel) and W = Dy/3/2 (lower panel). 
The density profiles (p*(y) = p(y)D) along y axis are shown in the insets. The 
inflection points are marked by diamond symbols. The pressure components are 
made dimensionless as follows P* t = (3PuDW, while the packing fraction is defined 
as 7] = py/D 2 - W 2 . 



Structural properties of hard disks in a narrow tube 



12 





Figure 4: Vertical correlation length (£* = £/D) as a function of linear pressure 
(P* = /3P XX WD) for W = L>/ 2 (upper panel) and W = 1)^/2 (lower panel). Vertical 
pair correlation function (g(i,j) = (yiUj)) are shown as a function of pair separation 
(\i — j\) in the insets. 




Figure 5: Pair distribution function (g(x)) as a function of distance (x* = x/D) for 
W = D/2 (upper panel) and W = Dy/3/2 (1 ower panel). MC simulation snapshots are 
shown at the corresponding packing fractions. 
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Figure 6: Equation of state (longitudinal pressure component vs. packing fraction) 
of the confined hard disks and that of analytical model where P* = (3P XX WD and 
r) = pVD 2 - W 2 . 

systems, the second peak completely disappears, the distribution function is extremely 
peaked and the distance between the peaks is less than D. This can only happen if the 
particles are in the alternating solid-like structure. These results are confirmed by the 
corresponding simulation snapshots presented in figure [51 too. Note that the locations 
of the inflection points in the transverse pressure are in accordance with the structural 
changes of the pair distribution function and the simulations snapshots. 

In the rest of this section our analytical and the numerical results are compared and 
we discuss the strength and the weakness of the applied approximation for the contact 
distance. Figure [6] shows the analytical and the numerical longitudinal pressures as 
a function of packing fraction for W = 0.8D. We can see that the low and high 
density behaviours of both solutions agree very well. This is not surprising because 
both are ideal gas at low density and both are packed in the same alternating (triangular 
lattice-like) structure at the close packing. However, the two pressure curves deviate 
substantially at the intermediate pressures and even a crossover takes place between 
them at 7] m 0.5. It is easy to understand that these behaviours stem from the 
exaggeration of the contact distance for small transverse positional differences. Figure 
[2] shows that the contact distance is overestimated by about 70% for \y,i — = 

in the case of W = 0.8D. Consequently the particles of the analytical model always 
interact with larger contact distances in all horizontal positions and they occupy larger 
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Figure 7: Comparison of the correlation lengths of the confined hard disks and that of 
the analytical model, where £* = C./D and P* = f3P xx WD. Solid curves corresponds to 
the analytical model, dashed curve to the asymptotic formula for the analytical model 
and short dashed curve for the hard disks. 



regions of the space. As a result, the analytical model possess higher pressure in the 
Tonks gas regime and it undergoes the fluid-solid structural change at lower densities 
than the system of hard disks. We can also see that the particles are pushed more 
strongly to the walls in the analytical model than in the numerical one to produce lower 
pressures at high densities. In the analytical model it is feasible to determine the high 
density limit of the pressures analytically. We get that /3P XX W = 2p/(l - pVD 2 - W 2 ) 

and P yy /P xx = |(l + 1/(1 — W 2 /D 2 )). The first expression tells us that the structural 
change gives an extra factor two in the well-known Tonks equations, while the second 
one shows that the transverse pressure becomes higher than the longitudinal one for 
all possible widths at the close packing. Note that the above relation between P xx and 
P yy is probably not true for hard disks and the phase behaviour of the two systems is 
not the same at high densities. We cannot prove it, because the eigenvalue equation 
(jSJ) cannot be solved numerically at very high pressures. However the examination of 
the correlation length reveals the differences residing in the two models (see figure EJ. 
For the analytical model an expression can be derived for very high pressures, which 
is given by £ = — exp(— c(3P xx W 2 )/(4c(3P xx W 2 ), while £ is determined from equations 
(1 1 5 a[) and fU6al) . Apart from the low-density limit the numerical and the analytical 
correlation lengths diverges with different slope with increasing pressure. We can see 
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Figure 8: Density profile of confined hard disks (solid curve) and that of analytical model 
(dashed curve) at pore width W = D/2 for r) = v / 3/2. Inset shows the corresponding 
transverse pair correlation function, where diamonds correspond to the confined hard 
disks, while the filled circles to the analytical model. 

that the analytical model is more correlated than the numerical one at the same pressure. 
It is also clear that the pressure is sufficiently high for the comparison because the 
analytical expression for the high pressure limit of the correlation length fits the exact 
one. The reason for the difference in the high pressure limit of the correlations is that 
the analytical model can pack more efficiently than the numerical one. The numerical 
model corresponds to the confined system of hard rhombuses with longer horizontal 
extension, i.e. it is harder for them to move between the two walls without overlapping 
with the neighbors. The uncorrelated movement is much easier for the disks, because 
they are more rounded than the rhombuses. Therefore the first order expansion of 
the contact distance gives rise to a stronger tendency for positional ordering and the 
two systems do not behave identically for very high pressures. Finally we compare the 
density profiles and the correlation functions of the two models in figure [SJ It confirms 
our present results that the analytical model is more ordered than the hard disks at 
the walls, the correlations are also stronger in the analytical model and the structure of 
solid-like clusters are triangular lattice-like in both models. 

In summary we can say that apart from some quantitative differences between the 
two models, it is proved very useful to devise the analytical model because it follows 
closely the phase behaviour of the confined hard disks. Furthermore analytical equations 
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can be derived for the pressure, density profile and correlation function, which make 
possible to examine the high pressure limit of the equation of state and the correlation 
length. 

5. Conclusions 

In this paper we have studied the phase behaviour of two-dimensional hard disks confined 
between two parallel lines by the transfer-matrix method and MC simulations. The 
width of the tube is restricted to the range of < W < Dy/3/2 to guarantee the 
exact statistical mechanical treatment of the system by taking into account the nearest 
neighbour interactions only. In summary the following important results are obtained: 

(i) The numerical investigation of the thermodynamic and structural properties has 
shown that no phase transition occurs in very narrow tubes. However there are 
some indications of the structural change between a fluid-like and the solid-like 
structures. These are the followings: 

(a) the longitudinal pressure component may possess two inflection points, 

(b) the longitudinal pressure deviates from the ideal gas law and it may exceed 
the transverse pressure component, 

(c) the pair correlation function shows damped oscillation, decays exponentially 
with the distance and it becomes infinitely long-ranged in the close packing 
limit, 

(d) the pair distribution function changes from the typical Tonks gas behaviour 
to a solid-like one having shorter period. 

(ii) It has been managed to device an exactly solvable system, which is obtained 
by the first order Taylor expansion of the contact distance between two hard 
disks. It has been shown that the new system corresponds to the system of 
parallel rhombuses between parallel walls. The eigenfunctions and eigenvalues and 
the corresponding pressure components, density distribution and pair correlation 
function are determined in analytical form. The analytically solvable model 
reproduces all important properties of the confined hard disks and makes possible 
to determine the close packing properties of the pressure and the correlation length. 

(iii) The numerical, analytical and simulation studies show together that the disks 
arrange themselves into an up-down sequence, where the neighboring particles 
are forced to stay at the opposite walls. The size of the alternating (triangular 
lattice-like) clusters becomes larger with increasing density and diverges at the 
close packing density. At such packing where the dominant interactions changes 
their directions from the horizontal direction to the vertical one, the longitudinal 
pressure and the pair distribution function shows unexpected behaviours. 

Even though our present study is exact only for very narrow tubes we believe that it is 
a step forward to understand the crystallization of two-dimensional bulk and confined 
fluids. It would be especially interesting to go beyond the nearest neighbor restriction 
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and examine the effect of the width on the nature of the structural rearrangement. It 
may happen that the correlation length diverges at finite density and first order phase 
transition takes place between a fluid and solid phases. In very narrow pores it may 
possible to examine the dynamical properties of the hard disk analytically, which may 
give deeper insight into the triangular lattice solidification of the fluid. One step ahead in 
this direction is the recent study of Arenzon et. al. |34j . where the dynamical properties 
of the one-dimensional system of rotating hard rectangles have been examined. In the 
very recent study of Herrera- Velarde et al. [35], common structural and dynamical 
behaviors have been found at the transitions from fluid-like to solid-like structures in 
some repulsive one-dimensional systems. 
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